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Abstract 

We present an improved version of the dual parton model which contains 
a new realization of the diquark breaking mechanism of baryon stopping. We 
reproduce in this way the net baryon yield in nuclear collisions. The model, 
which also considers strings originating from diquark-antidiquark pairs in the 
nucleon sea, reproduces the observed yields of p and A and their antiparticles 
and underestimates cascades by less than 50 %. However, ri's are under- 
estimated by a factor five. Agreement with data is restored by final state 
interaction, with an averaged cross-section as small as o" = 0.14 mb. Hyperon 
yields increase significantly faster than antihyperons, in agreement with ex- 
periment. 
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1 Introduction 



A striking feature of heavy ion collisions is the huge stopping of the partici- 
pating nucleons. At CERN energies, the deep minimum in the net baryon rapidity 
distribution {AB = B — B) at y* ^ 0, observed in pp collisions, has been practically 
filled up in a central collision of heavy ions ^ . For central Pb-Pb collisions, the 
value of this density is five times larger than the corresponding value in pp - scaled 
by the average number of participants. Note that the total number oi B — B (i.e. 
integrated over rapidity), exactly satisfies scaling in the number of participants, due 
to baryon number conservation. This shows the dramatic change in the shape of the 
B — B (and p — p) rapidity distributions between pp and central Pb-Pb collisions. 
Such a change is usually referred to as baryon stopping. 

All independent string models of hadronic and nuclear collisions in their original 
form completely fail to reproduce this important feature of heavy ion collisions. In 
the dual parton model (DPM) and in the quark gluon string model (QGSM) |Q, 
the dominant contribution to particle production in pp collisions at ~ 20 GeV, 
consists in two qq-q strings, which produce, after fragmentation, two baryons in the 
fragmentation regions of the colliding protons. Starting with the Lund model, which 
initially had a single string, the above mechanism of particle production has been 
adopted in most current string models. In these models, there is some amount of 
stopping due to energy conservation. This produces an increase of the net baryon 
yield at mid rapidities between NN and central Pb-Pb collisions, which is typically 
of a factor 2 - more than two times smaller than the observed one. Hence, the 
dramatic failure of all these models to reproduce the observed stopping. 

Actually, a possibility to slow down the net baryon in pp collisions, was intro- 
duced a long time ago by Rossi and Veneziano 0. In their approach, the baryon is 
viewed as three valence quarks bound together by three strings each one with a quark 
at one end and with the other end joining in a point called string junction. This 
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string junction carries momentum as well as the baryon quantum number. Rossi 
and Veneziano pointed out that the string junction could migrate to mid rapidities 
with a distribution in da/dx ~ l/v^ (or da/dy ~ exp(— l/2|(y — ymax\)- This cor- 
responds to an annihilation cross-section which decreases with energy like s'^^"^. In 
refs. [0 and a distribution of the string junction in (i.e. flat in rapidity) was 
proposed, corresponding to an annihilation cross-section which reaches a constant 
asymptotic value (of 1 to 2 mb). Here we adopt the first approach. However, we do 
not rule out the second possibility - which would have important consequences at 
the energies of the future heavy ion colliders p. 

The above stopping mechanism has been recently introduced in heavy ion col- 



lisions ^ and implemented in the Hijing and Venus [|Tl|] Monte Carlos. 
However, the introduction of the Rossi- Veneziano mechanism does not explain by 
itself why the stopping is larger in central heavy ion collisions than in pp. In ref. 
1^, a mechanism to enhance stopping in heavy ion collisions was proposed. It was 
based on the separation of the pp cross-section app = a^J' + into a diquark 
breaking {DB) and a diquark preserving {DP) piece, and on the assumption that 
the diquark can be broken in any inelastic collision. These results in a DB cross 
section in pA and AA collisions which increases faster with A than the DP one. The 
drawback of this approach is that it requires some fine tuning. The value of has 
to be small enough in order not to contradict the pp and pA data (where stopping 
is comparatively small) and large enough to produce the large stopping observed in 
central heavy ion collisions. 

In a recent publication |12|, a new formulation of the DB mechanism has been 
introduced in which this drawback is avoided (i.e. one can have a^^ negligeably 
small at CERN energies, and still have an important effect in central Pb-Pb colli- 



sions). In the present paper we use the formulation of |JT2[ to compute the rapidity 
distributions of i? — 5 in hadronic and nuclear collisions. We obtain a reasonable 
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agreement with experiment0. 

Another striking feature of the CERN heavy ion program is the strong increase 
of the yields of hyperons and anti-hyperons per participant between pp or pA and 
central AB collisions. This increase obeys to the hierarchy Q > S > A (i.e. the 
large number of strange quarks in Y, the larger the increase) [|T3], In two 

recent publications |]T3[, it has been shown that a baryon stopping mechanism 
of the type described above produces a substantial increase of the hyperon yields 



according to this hierarchy. The physical reason for this increase is quite obvious |16 



Since in the diquark breaking component the net baryon is formed out of three sea 
quarks around the string junction, the probability of producing hyperons is strongly 
enhanced - especially for f2's since its probability of production in the conventional 
diquark fragmentation mechanism is zero. In the present paper we extend the results 



of 1 12] in two directions. First, we study the rapidity distributions of protons and 



hyperons in pA and AA collisions (in refs. [0 and [15| the analysis was restricted 
to the rapidity window \y*\ < 0.5). Second, we show how the four free parameters 
of ref. [|T^ can actually be reduced to two. This makes the model more predictive, 
especially for the antihyperon over hyperon ratios. 

While the yield of A's and, to a large extent, of cascades can be described by 
the model, that of omegas is underestimated by almost an order of magnitude. The 
same conclusion has been reached in [O, |15|1. We show that final state interaction, 



with an averaged cross-section as small as cr = 0.14 mb allows to describe all 
hyperon and antihyperon yields. 



A similar value of a was found in ref. in the hadron gas model. It was 
argued there that, due to this small value of a, interactions in a hadron gas could 
not drive the system to chemical equilibrium (the process would be too slow). We 
find, indeed, that the effect of final state interaction in p and A production is very 



''Note, however, that similar results are obtained using the approach of ref. g| with the DB 
mechanism of Fig. 2. 
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small. Its effect on S production is moderate. Only for such a rare process as 
production its effect is very important - making the f2 + yield five times larger 
than the value obtained without final state interaction. 

Due to the diquark breaking component, we obtain an increase of hyperons 
substantially larger than the one of antihyperons, i.e. the ratio between Pb-Pb and 
pPb yields is substantially larger for Y than for Y. This effect is enhanced by final 
state interaction. It has been observed experimentally jl^ . 

The paper is organized as follows. In Section 2 we describe the baryon stopping 
mechanism and compute the net baryon {B — B) rapidity distributions in pp, pPb, 
SS and Pb-Pb collisions. In Section 3 we compute the rapidity distributions of p— p, 
and Y — Y . In Section 4 we describe BB pair production from strings containing 
sea diquarks or antidiquarks at one of their ends and show how the y4-dependence 
of B production is increased. In Section 5 we study the effect of the final state 
interaction, separately on the Y and Y yields. Section 6 contains a discussion of our 
results. Conclusions are given in Section 7. 

2 Baryon stopping 

A fragmentation string mechanism in which the x — > and x — > 1 behaviour of 
the fragmentation functions is controlled by Regge intercepts has been introduced 
in ref. [0. In the case of net baryon production, it consists of a sum of two 
terms as depicted in Fig. 1. In the original Lund fragmentation scheme [|l9l only 
the first term is considered. The second one was introduced later - the so-called 
pop-corn mechanism. Even with the inclusion of the second component (Fig. lb), 
this mechanism leads to the production of too fast baryons and fails completely to 
reproduce the observed stopping in heavy ion collisions. This fragmentation scheme 
(including the component in Fig. lb) will be referred to as the conventional or 



diquark preserving {DP) mechanism. Following ref. [O] we introduce the baryon 
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stopping mechanism showed in Fig. 2. It will be referred to as the diquark breaking 
(DB) component. In this component, the rapidity distribution of the produced net 
baryon AB = B — B in a N-N collision is 

(1 - + zi/' (1 - z_)"'-''/^] (1) 

where Z± = exp{±y — ymax) , n-i and n2 are the average number of collisions suffered 
by the two colliding nucleons and C„j^„2 is determined from the normalization to 
two. The factor Z^/^ has already been discussed in the Introduction. The factor 
{1 — Z) gives the behavior near y = y-max- There is some uncertainty concerning its 
power 1^. The value in eq. (|I|) is obtained as follows. From Fig. 2 we see that in 
order to produce the baryon at y ~ ymax it is necessary to slow down three quarks. 
Assuming they behave as 1/y/x at the energies under consideration 0], we obtain a 
power 1/2 for the case of Fig. 2 - which corresponds to n = 2. In the general case 
of n inelastic collisions we obtain the power n — 3/2 in eq. (0). 
The corresponding distribution in AA collisions is fT^ 
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Here ua and n are the average number of participants of nucleus A and the average 
number of collisions, respectively. dNijp/dy is given by the conventional, diquark 
preserving, hadronization mechanism for which we use the results of ref. [^, and 
dNoB/dy is given by eq. ([I|). The integral over y of both rapidity distributions is 
equal to two (baryon number conservation). 

Let us discuss the physical meaning of eq. @ in pp interactions (n^ = 1). We 
see that in the case of a single inelastic collision {n = 1), we recover the conventional 
DP mechanism. The underlying assumption is that, in this case, the string junction 
follows the valence diquark and baryon production takes place in the conventional 
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way. (We do not exclude a small admixture of the DB component in this case, 
but experimental data do not require its presence). Consider next the case of two 
inelastic collisions {n = 2). Here the underlying assumption is that there is an equal 
probability (1/2) for the net baryon to be produced in any of the two collisions. 
However, in only one of them can the string junction follow the valence diquark and 
fragment in the conventional {DP) way. In the other one, the string junction is 
free and baryon production takes place according to the DB mechanism. The latter 
is responsible of (most of) the observed baryon stopping. The generalization to n 
inelastic collisions and to A A interactions (eq. @)) is then straightforward. 

It is quite remarkable that such a simple mechanism, with no free parameter 
(modulo the uncertainty in the power of 1 — Z in eq. (|lD, discussed above) gives a 
good description of the present data on the net baryon rapidity distribution. 

Note that in the case of pp interactions at i/i ~ 20 GeV, the two string compo- 
nent [n = 1) dominates, and, as discussed above, we recover the usual DP results. 
With increasing energies, the components with n ^ 1 become increasingly important 
and baryon stopping will increase. Eq. (|I|) (with = 1) not only gives definite pre- 
dictions concerning this increase but, moreover, leads to specific qualitative features. 
In particular stopping will strongly depend on the charged particle multiplicity. A 
low multiplicity event sample selects low values of fi where stopping will be compar- 
atively small, while, at large multiplicities, stopping will be larger. Such a feature 



has been observed recently at HERA - and discussed in ref. |22l] in a different 
theoretical framework. 

We turn next to the generalization of eq. (H) to asymmetric interactions such as 
pA. In this case one has 
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Here dN'^'^^~'^^^'^^^ /dy denotes the rapidity distribution of a string stretched between 
a diquark of one of the n wounded nucleons of A and a valence (or sea) quark of 
the proton. (It is computed in DPM as a convolution of momentum distribution 
functions and fragmentation functions). Since each of the wounded nucleons suffers 
a single inelastic collision, only the DP component is involved - with each diquark 
fragmenting in the nucleus fragmentation region {y* < 0). The terms in the bracket 
correspond to the fragmentation of the incoming proton. Since it suffers n inelastic 
collisions, we have the DP hadronization mechanism (with probability 1/n) and the 
DB one (with probability (n — l)/n). The latter is now given by the first term of 
eq. (0). All rapidity distributions, integrated over y, are equal to one in this case. 

The B ~ B rapidity distributions obtained from eq. (^ in central SS and Pb-Pb 
collisions are shown in Fig. 3. Note that at mid-rapidities, these distributions are 
dominated by the DB component. Not only the latter is proportional to n — fiA (eq. 
(ID), but, moreover, dN^s/dy is larger than dN^p/dy at mid rapidities (see Table 
1). Nevertheless, the existence of the two huge maxima of the DP component in the 
fragmentation regions, still shows up in the AA distribution. For a given system, the 
detailed shape of the B — B rapidity distribution depends on the power of {1 — Z) 
in eq. (p. As discussed above there is some theoretical uncertainty in the value of 
this power. However, the variation of the shape of this rapidity distribution from 
one system to another is a characteristic feature of the model. As seen in Fig. 3, the 
minimum at mid-rapidities is gradually filled up from pp to central Pb-Pb collisions 
and, therefore, it is more pronounced in SS than in Pb-Pb. 
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3 Net hyperon enhancement 

In the previous section we have shown that baryon stopping can be described 
using a new formulation of the diquark breaking (DB) mechanism. In this case, 
depicted in Fig. 2, the string junction is surrounded by three sea quarks to produce 
the net baryon. Therefore, not only the net proton yield Ap = p~p will be strongly 
enhanced from pp to central AA collisions, but also the net hyperon yield Ay = 
Y ~Y. This is especially so for fi's, which cannot be produced at all with the DP 
mechanism of Figs, la and lb[]. More precisely the ratio between yields in central AA 
and pA (or pp) collisions will obey to the hierarchy : AQ > AH > AA > Ap. This 
is in agreement with the results of the WA97 [14] and NA49 [|l], |13| collaborations. 

The relative yields of the different baryon species will be determined by the 
strangeness suppression factor S/L where S is the probability associated to the 
strange quark and L the one associated to the light quarks {u or d). We consider 
two possibilities : S = 0.10 and L = (1 - S)/2 = 0.45 {S/L = 0.22) and S = 0.13 
and L = 0.435 {S/L = 0.3). With baryons produced out of three sea quarks (Fig. 2) 
it is easy to see that the relative yields are 

h = 4:L^ : 4L3 : 12^^^ : 3LS^ : 3LS^ : S^ (4) 

for p, n, A + S, and Q, respectively. Moreover, we take, S+ + = 0.6A. 

This reduction in the number of charged S's is due to resonance decay (S(1385)Pi3 
decays into Ait with an 88 ± 2 % fraction). 

It is interesting that, in spite of the huge hyperon enhancement observed experi- 
mentally, the factors (^ lead (both with S = 0.10 and S = 0.13) to an overestimation 
of hyperon production in pPb collisions especially for H's and Q^s. In central Pb-Pb 



collisions, net hyperon production is also overestimated - except for f2's. In ref. 112 1 



new component consisting in a diquark which contains sea quarks has been introduced in 
p3|] . However, this diquark is assumed to have the same momentum distribution as a diquark 
made out of two valence quarks and, hence, produces baryons mainly in the fragmentation regions. 
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the following explanation of this hyperon excess was proposed: at present energies, 
it may happen that the net baryon is not formed out of three sea quarks as in Fig. 2, 
but, due to phase space limitation, a valence quark, at one of the ends of the string 
where the baryon is produced, is picked up together with two sea quarks. Obviously, 
in this case the strangeness production rate is substantially reduced (in particular, 
Q production is not possible in this case). The relative yields J3 in (^) are then 
changed into 



In the following, we introduce a free parameter a (0 < a < 1) which determines the 
admixture of and I2 given by (^) and (|^). More precisely, we will take the relative 
yields given by 



The best description of the data is obtained with a = 0.23 for 5 = 0.13 and a = 0.5 
for S = 0.1. The results are similar in the two cases ; the sensitivity to the value 
of the strangeness suppression factor S/L turns out to be quite small. The results 
for the rapidity distribution of the net yields of p and A in central Pb-Pb collisions 
are given in Figs. 4 and 5. Note that our rapidity distribution for A — A (Fig. 5) is 
broader than the estimates of the NA49 collaboration |l[]. This, in turn, produces 
some discrepancies in the p — p yield (Fig. 4). Final data on A and A are needed in 
order to clarify the situation ; we shall come back to this point in section 6. The net 
A rapidity distribution in central SS collisions is shown in Fig. 6. The corresponding 
results for minimum bias pA collisions are given in Figs. 7 and 8. It is seen that the 
normalization of the experimental data is larger than the theoretical one especially 
for p — p. Note, however, that by integrating over y the experimental distributions 
one realizes that their normalization is larger than the number of wounded nucleus 



I2 = 2L^ ■ 2L^ ■ ALS : : 8^2 : 



(5) 



/ = a J3 + (a - 1)J2 



(6) 
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in pAu (given by the Glauber model) by more than a factor 2. As pointed out in 
this excess may be due to recoil nucleons which are not completely disentangled 
from the wounded ones. This point needs clarification. 

Note that, in our model, the relative yields I3 (eq. (^)) should apply at higher 
energies when the phase space limitations are less important. With S/L = 0.3, they 
would give a ratio H+ + E'^/A + S ~ 0.3 which is in agreement with Fermilab 
and SPS collider data 



4 Antibaryon production 

In string models, BB pair production takes place via diquark-antidiquark pair 
production in the string fragmentation. It turns out that at present CERN energies 
only strings of type qq-q have large enough invariant mass to produce BB pairs. This 
gives rise to a scaling of B yields in the number of participants. (The number of qq-q 
strings is proportional to the number of wounded nucleons). Experimentally, the 
observed increase is much faster - closer to a scaling in the number of collisions. In 



order to solve this problem it was proposed some time ago |^0|, ^ to consider the 
production of BB pairs from diquark-antidiquark pairs in the sea of the participating 
nucleons[|. The rapidity distribution of antibaryons in AA collisions is then given by 

^ {y) = nA\—^ — -{y)\ +{n-nA)[—^{y)] ■ (7) 

The first term represents the conventional pair production in the string breaking 
process. As discussed above, it is proportional to the number of participants. The 
second term corresponds to pair production from a string having a sea diquark or 
antidiquark at one of their ends. In DPM, the total number of strings is proportional 
to n. Since the number of strings with a valence diquark at one of their ends is 



"a different mechanism based on string junction-antijunction exchange has been proposed re- 
cently 
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proportional to ha, the number of strings with a sea diquark at one of their ends is 
proportional to n—fiA- Of course, pulling a diquark-antidiquark pair out the nucleon 
sea is dynamically suppressed - in the same way as its production in the string 
breaking process is suppressed as compared to q-q production. Thus, we expect 
that in each individual string, the production of BB pairs in the two components 
(sea and string) in eq. (J^) are comparable. In practice, the normalization of the 
second component is treated as a free parameter. However, it turns out that the sea 
component is always smaller than the string one, not only at = (see Table 1) 
but also after integration over rapidity. Note that the string with sea diquarks have 
a smaller invariant mass. 



For the string term we use the results of ref. The absolute normalization 

of this term was determined from a fit of the pp data. For the y-dependence of 



the sea term, we also use the results of ref. [0. As discussed above, its absolute 
normalization is a free parameter. This parameter is the same for all species of 
baryons. More precisely, since the baryons and antibaryons in the sea component 
are made out of three sea quarks or antiquarks, the relative yields of the different 
baryon species is again given by eq. (^)[^. (We neglect here the small differences in 
the rapidity shapes induced by the different baryon masses) . We are left in this way 
with a single free parameter for this new sea component. Therefore, we have a total 
number of two free parameters, one in the diquark breaking component and one in 
the sea component - plus the value of the strangeness suppression factor S/L for 
which two values (0.22 and 0.3) have been considered. Of course, the conventional 
components DP and string in eqs. (0) and (^ contain several free parameters. 



However, as discussed above, these parameters have been fixed in ref. pO| from a 

fit of the pp data and are not changed here f^. The values aX y* = of the various 

**Note that in this case a — 1, i.e. no admixture of the type discussed in connection with the 
DB component is present here. 

^'^For this reason, the change in the strange suppression parameter S/L from 0.22 to 0.3 only 
apphes to the new components DB and sea in eqs. (||) and (0) (see Table 1). 
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components for the different baryon species are given in Table 1. 

The resuhs for the rapidity distributions of the A + A, S~ + S~ and + H 
yields in central Pb-Pb collisions at 158 GeV are given in Figs. 9-12. The p, p, 
Y, and Y yields at \y*\ < 0.5 in minimum bias pPb collisions as well as at four 
different centralities in Pb-Pb collisions, are given in Fig. 13. The corresponding 
ratios Ry = Y/Y at \y*\ < 0.5 are given in Fig. 14. It should be noted that 
the values of Ry are not absolute predictions of our model. They can be changed 
by changing the normalization of the sea component in eq. (0). However, the 
ratios Rp : R\ : R^ : Rq are a characteristic feature of the model. They show an 
increase with the number of strange quarks in the baryon. In Fig. 15 we show the 
various baryon yields at y* = divided to the average number of participants, 2n^, 
normalized to the same quantity in pPb. A discussion of these results is given in 
Section 6, after introducing final state interaction. 

5 Final state interaction 

In an attempt to explain the strong enhancement of the 1] + yield observed by 
the WA97 collaboration [|1^], we are going to use our results for the baryon densities 
as initial conditions in the gain and loss differential equations which govern final 



state interactions [17, 28 



^ = pk{x) pe{x) - J2 ^ik Piix) pk{x) . (8) 

The first term in the r.h.s. of (§) describes the production of particles of type 
i resulting from the interaction of particles k and £ with space-time densities p{x) 
and cross-sections a^i (averaged over the momentum distribution of the interacting 
particles). The second term describes the loss of particles of type i due to its 
interaction with particles of type k. We use cylindrical space-time variables and 
assume boost invariance (i.e. the densities p(x) are taken to be independent of y). 
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If we furtliermore assume tliat tlie dilution in time of the densities is mainly due to 
longitudinal motion, i.e. : 

pi{^,y,s) = pi{T,s) — , (9) 

r 

where r = y/t"^ — is the proper time and s the transverse coordinate, eqs. ® can 
be written as 




— * — # —^ 

Here Pi{y, s, h) = dNi/dydsdb. Thus, at fixed impact parameter 6, we have to know 
the rapidity densities per unit of transverse area ds. Our eqs. (0) and do give 
these rapidity densities - the dependence on s and b is contained in the geometrical 
factors fiA and n, given by the Glauber model. In the following, we use nuclear 
profiles obtained from Saxon- Woods nuclear densities using the three-parameter 
Fermi distribution of ref . . For the pion densities we use the DPM results of ref . 
pod - where explicit expressions as a function of and n are given. 

Eqs. ( |T0|) have to be integrated from initial time tq to freeze-out time r. These 
equations are invariant under the change r — > cr. Therefore the result depends 
only on the ratio t/tq. Following refs. 13^, we use the (inverse) proportionality 
between r and p and put t/tq = p{y, s, b)/ pf^y Here p(?/, s, b) are the initial densities 
given by our expressions obtained in previous sections and pf^ is the freeze-out 
density. For the latter, we take the charged density per unit rapidity in a pp collision, 
i.e. pf, = {3/7rRl]idN~/dy)y,=o = 1.15 fm'^ |30|, 0]. 

We have now to specify the channels that have been taken into account in our 
calculation. They are : 



TiN ^KA , ttN ^ , nA^KE , ^ KE , nE ^ KQ , (11) 



14 



and the corresponding reactions for antiparticles. To be more precise, of all possible 
charge combinations in (pT]), some are of the type shown in Fig. 16a, with anni- 
hilation of a light quark pair and production of an s-s. They have all been taken 
into account with the same cross-section a = 0.14 mb. All other reactions in ( ]TT| ) 
are of the type shown in Fig. 16b. They have three quark lines in the t-channel 
(baryon exchange). Their average cross-section are smaller than the one of Fig. 16a 
and have been neglected^. We have also neglected all strangeness exchange reac- 
tions (Fig. 16c) KN < — y ttA etc. Although the corresponding cross-sections are 
larger at threshold, this is no longer the case for the cross-sections averaged over 
the momentum distributions of the interacting particles . (This is due to their steep 
decrease with increasing energy (see |T^). Channels (pUj) are thus dominant due 
to the relations Pn > Pa > Ps > Pn and p,r > Pk between particle densities. The 
results, obtained after solving numerically eqs. (^), with our initial densities and 
a common value of the averaged cross-section a = 0.14 mb for all channels, are 
shown in our figures by a full line in the case S/L = 0.3 and a = 0.23 and by a 
dashed-dotted line in the case S/L = 0.22 and a = 0.5. (A comparable value of a 
has been obtained in ref. [|T^ in a hadron gas model). 

The effect of the final state interaction is negligeably small in pA collisions. In 
central 5*5* collisions its effect on the p and A yields is very small (less than 5 %). 
The effect increases with the number of strange quarks in the produced hyperon. 
In central Pb-Pb collisions, with our value of the cross-section, it turns out to be 
comparatively small for p and A yields. However, it increases the S yields by up 
to 50 % and the f2 + H yield by a factor 5. Agreement with the WA97 data fl^ is 

obtained in this way (Fig. 13). 

WThe reactions we have kept are : 7r+ + n ^ i^+A, n^p K^K, tt^ + n ^ K'^J^^ , n^p 
K+i:+, TT-A ^ K"E-, TT+A ^ K+E°, Tr+S" ^ K+E' , Tr-^+ K°E°, tt-EP -> K^Q. and 
7r+S~ K^Q, for the reactions initiated by 7r+ or tt^. For aU of them, as well as for the 
corresponding ones with antiparticles, we take a = 0.14 mb. The reactions initiated by 7r° are 
either of the type of Fig. 16a or Fig. 16b, depending on whether the uu or dd component of the tt*^ 
is considered. For this reason all these reactions have been included with cross-section a/2. 
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It is important to note that, due to the small value of cr, the final state interaction 
has an important effect only on very rare processes such as Q production. It cannot 
drive the system into chemical equilibrium - even locally. 



6 Discussion 

We discuss here the main features of our results. For A + A our results for 
Pb-Pb are slightly higher than the WA97 data and grossely underestimate NA49 
ones at mid-rapidities (Fig. 9). Note, however, that the latter are very preliminary 
and are currently under reanalysis. For central SS collisions, where the NA35 data 
are final, we slightly underestimate their net A yield and slightly overestimate the 
total A yield form NA36 (Fig. 6). However, the NA35 value for the A yield at 
mid rapidities 0.75 ± 0.15 is about two times larger than our result. Note that this 
experimental point looks "anomalous" : compared with the WA97 value for the most 
central rapidity bin in Pb-Pb collisions, (1.8 ± 0.2), there is an increase by a factor 
2.4 - whereas the number of participants increases by a factor 7. Note also that 
NA35 finds a ratio A/p = 1.9±0.7at mid-rapidities, while in our model this ratio 
is always smaller than one (see Fig. 13). This important point needs clarification. 
In particular, final values of this ratio in Pb-Pb collisions are needed. 

Concerning, the cascade yields, the NA49 data are published fl^. They are 



30 40 % higher than the WA97 ones ||Tj] at mid rapidities. Our results, after final 



state interaction, are intermediate between the two sets of data - but somewhat 
closer to the NA49 results (Figs. 10 and 11). The 17 + yields are in agreement 
with the WA97 data, after final state interaction (Fig. 12 and 13). 

As discussed in Section 4, in our model, the ratios Ry = Y/Y increase with the 
number of strange quarks in the baryon (Fig. 14). This tendency is also seen in the 
data. However, the ratio of ratios R-^/Rk is somewhat too small in our model as 
compared to the WA97 data |[T^, but agrees with the NA49 ones [|T], [13]. (Remember, 
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however, that the value of i?A from NA49 is prehminary) . 

Another characteristic feature of our approach is that, at mid-rapidities, hype- 
rons are more strongly enhanced than antihyperons. As a consequence, the ratio 
Y /Y decreases between pPh and central Ph-Ph collisions (Fig. 14). This is due to the 
strong effect of the DB component in the net baryon yield. Final state interaction 
works in the same direction. This important feature of our results is seen in the 
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Finally, the WA97 collaboration has found that the increase of the hyperon and 
antihyperon yields per participant (Fig. 15) increases faster than the number of 
participants between pPh and the first centrality bin in Ph-Ph. However, between 
the first and last centrality bin all yields approximately scale with the number of 
participants. We find an increase which is faster in the first case than in the second 
one. However, some mild increase is left in Ph-Ph (Fig. 15). 

7 Conclusions 

The large baryon stopping observed in central heavy ion collisions at CERN ener- 
gy is not reproduced by any of the available independent string models - at least in 
their original form. We have modified the DPM by introducing a new realization of 
the diquark breaking mechanism. We reproduce in this way the observed net baryon 
yield. This mechanism also produces an important enhancement of net hyperons. 
At this level, the new version of DPM presented here (which has also diquark- 
antiquark pairs in the nucleon sea) remains strictly an independent string model. 
It reproduces with two free parameters, the observed yields of p and A and their 
antiparticles in pA and Ph-Ph collisions. Cascades in central Ph-Ph collisions are 
underestimated by less than 50 % while f2's are too small by a factor 5. Agreement 
with experiment is restored by introducing final state interaction with an averaged 
cross-section as small as a = 0.14 mb. In this way, we depart from string indepen- 
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dence. However, with this small value of the cross-section, there is no significant 
efi^ect on the bulk of particle production. A comparable value of the cross-section 
for final state interaction was obtained in ref. Jl^ from the experimental data on 
the energy dependence of cross-sections, averaged over the momentum distribution 
of the interacting particles obtained in the hadron gas model. The smallness of this 
averaged cross-section lead the authors of [|T^ to argue that strangeness phase space 
saturation would be too slow in a hadron gas. It is interesting that such a small 
value of the averaged cross-section allows to reproduce the observed enhancement 
of multi-strange hyperons and antihyperons in central Pb-Pb collisions. 

The main features of our results are the following : 1) The hyperon yields per 
participant increase faster than antihyperon ones. As a consequence, the ratio 
Ry = Y/Y decreases between pPb and central Pb-Pb collisions ; 2) The ratios 
Ry increase with the number of strange quarks in the hyperon ; 3) The increase of 
the Y and Y yields per participant is faster between pPb and the first centrality 
bin in Pb-Pb collisions and slows down between the first and last centrality bins of 
WA97. All these features are also present in the data. 
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Figure Captions 



Fig. 1 : Conventional diquark preserving (DP) fragmentation mechanism for net 
baryon production. 

Fig. 2 : Example of diquark breaking (DB) diagram for net baryon production in 
pA with two inelastic collisions. 

Fig. 3 : Rapidity distribution of the net baryon number [B — B) in central SS 
(200 AGev/c) and PhPh (158 AGev/c) coUisions. The full hues are obtained 
from eq. (2). The data are from ref [1,2] (for central SS coUisions the data are 
obtained as B — B = 2 {p — p) + 1.6 (A — A)). Open circles arc data reflected 
about y* — (errors not shown) . The dotted line is the result obtained without 
the DB component in the case of central PbPb coUisions. Due to baryon 
number conservation, these results are not affected by final state interactions. 

Fig. 4 : Rapidity distributions for net proton production {p — p) in central PbPb 
collisions at 158 AGev/c compared to the results of ref. [1]. Open circles 
are data reflected about y* — (errors not shown). The dashed line is our 
result without final state interactions with a strangeness suppression factor 
S/L = 0.3, and the full line is the corresponding result with final state interac- 
tions. The dotted-dashed fine corresponds to a suppression factor S/L — 0.22 
and with final state interactions. The dotted line is our result without DB 
component and without final state interactions. 

Fig. 5 : Same as Fig. 4 for A — A. Now the three dotted lines are estimates by the 
NA49 Collaboration [1]. The experimental point a,t y* — is form the WA97 
Collaboration [14]. 
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Fig. 6 : Same as Fig. 5 for central SS collisions. The data are from NA35 [2] 
(circles). Also shown, the value of the total A yield measured by NA36 [32] 
(squares). The theoretical curves are computed with: a — 0.23, S/L — 0.3 
(full line); a = 1, S/L = 0.3 (dashed line); a = 0.23, S/L = 0.22 (dotted line) 
and a — 1, S/ L — 0.22 (dashed-dotted line) 

Fig. 7 : Same as Fig. 4 for minimum bias pAu interactions. The experimental 
data are from ref. [2]. 

Fig. 8 : Same as Fig. 5 for minimum bias pAu interactions. The experimental 
data are from ref. [2]. 

Fig. 9 : Same as Fig. 4 for A + A. The point at y* = (circle) is from the WA97 
Coll. [14]. The squares are from NA49 [1]. 

Fig. 10 : Same as Fig. 9 for + S+. The NA49 data are from ref [13]. 

Fig. 11 : Same as Fig. 9 for S". The NA49 data are from ref [13]. 

Fig. 12 : Same as Fig. 9 for + Q. The experimental point is from WA97 [14]. 

Fig. 13 : Yields of p. A, S", Q+Q, p, A and S+ for minimum hiaspPb (158 Gev/c) 
and central PbPb collisions (158 AGeV/c) in four centrality bins. Experimen- 
tal data are from WA97 [14] (black points) and NA49 [13] (open square). 
Full (dashed) hues are our results with (without) final state interactions for 
strangeness suppression factor S/L — 0.3. The dashed-dotted lines are our 
results with final state interactions for S/L = 0.22. 

Fig. 14 : Ratios B/B at y* = for minimum bias pPb and PbPb collisions 
in four different centrality bins at 158 AGev/c. Black circles, squares and 
triangle correspond to experimental data of WA97 [14] for A/A, S+/S~ and 
fl/fl respectively. Open squares are NA49 data [13]. Full lines are our results 
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with final state interactions and dashed hnes without final state interactions, 
both for S/L = 0.3. 

Fig. 15 : Same as Fig. 14 for the baryon yields at = in PbPb divided by 
the number of participant nucleons relative to the same ratio in minimum bias 
pPb. 

Fig. 16a : Quark diagrams for reactions (11) with light quark pair annihilation 
and s — s quark creation. 

Fig. 16b : Quark diagram for reactions (11) with three quark exchange in the 
t-channel 

Fig. 16c : Quark diagram for strangeness exchange reactions. 
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Table 





P 


A 




Q 


dNgJdy 
dN^E/dy 


4.80 X 10-3 
9.12 X 10-2 


2.15 X 10-3 
3.00 X 10-2 


3.21 X 10-^ 
2.86 X 10-3 


3.20 X 10-5 
1.23 X 10-^ 


dN^J/dy 

dN^tring/dy 


6.90 X 10-2 
8.50 X 10-3 


1.40 X 10-2 
2.26 X 10-3 




1.65 X 10-^ 




5.07 X 10-6 


dNfjdy 

dN^E/dy 


6.54 X 10-3 
9.56 X 10-2 


2.43 X 10-3 
2.62 X 10-2 


2.98 X 10-^ 
2.29 X 10-3 


2.4 X 10-5 
1.23 X 10-4 



Table 1 : Values of the rapidity densities at = in eqs. (0) and (|^ for central 
PbPb collisions (n^ = 178, n = 858) with a = 0.23 and S/L = 0.3 (first four lines) 
and with a = 0.5 and S/L = 0.22 (last two lines). The DP and string contributions 
are the same in both cases. The value of the DP component for 3" is not exactly 
0, due to the fragmentation mechanism of Fig. lb. However its value is very small 
as compared to the other components and has been neglected. 
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